Inducing topological order in dirty wires: Majorana fermions from scattering 
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We focus on inducing topological state from regular, or irregular scattering in (i) p-wave super- 
conducting and (ii) proximity coupled Rashba wires. We find that while disorder is detrimental 
to topological state in p-wave wires, we find that it can induce topological state in Rashba wires 
contrary to common expectations. We find that the total phase space area of the topological state 
is conserved for long disordered wires, and can be even increased in an appropriately engineered 
superlattice potential. 
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The recent observation of a zero-bias peak in the An- 
dreev conductance of superconducting InSb nanowire 
heterostructures [l] has been attributed to Majorana 
fermions arising in quantum wires with effective topo- 
logical superconducting order [2^4] . Similar observations 
were consequentially reported by other groups |6]. 
Apart from their fundamental importance (being parti- 
cles that are their own antiparticles), Majorana fermions 
are also of technological importance as they provide 
a physical setting for topological quantum computing 
which is expected to be resistant to decoherence [2^|71|8]. 

In condensed matter settings Majorana fermions ap- 
pear in the form of zero-energy Andreev bound states in 
the so-called topological superconductors. Such topolog- 
ical superconducting order can arise in spin-orbit cou- 
pled semiconductor quantum wires in proximity to con- 
ventional superconductors [31 14] . Following the original 
proposals, which assumed single channel, disorder free 
nanowires, there was a rush of activity on whether the 
topological state, and thus the Majorana fermion would 
survive in multichannel wires jOHTT] or the presence of 
disorder p!2HT6] . The main conclusion of these works 
is that the topological state and hence the Majorana 
fermion survive provided (i) the mobility is high enough 
such that the localization length is shorter than the co- 
herence length of the topological superconductor and 
(ii) there is an odd number of spin-resolved transverse 
modes in a multi-mode wire. A side conclusion of (i) is 
that the disorder is always detrimental to the topological 
order. 

In contrast, the present work suggests that both con- 
ditions (i) and (ii) are in fact not necessary for having 
Majorana fermions. In particular, we show below that 
disorder can induce topological order and thus create 
Majorana fermions, or even facilitate their experimen- 
tal realization. The apparent contradiction with regard 
to previous work is resolved by the fact that many of 
the previous results were based on the effective model 
of a spinless p-wave Hamiltonian (arising in the contin- 
uum limit of he Kitaev chain [2]). This model arises as 
an approximate Hamiltonian for the experimentally rel- 



evant proximity-coupled semiconducting nanowire in the 
almost depleted limit. Below, we go beyond this descrip- 
tion to a model that treats spin-split bands realistically 
and show how random (disorder) or regular (e.g. due to 
a superlattice) scattering can create topological order. 

In this Letter, we therefore discuss the effects of regular 
and irregular scattering on the topological state. To this 
end, we develop a theory capable of studying topologi- 
cal phase transitions in the presence of individual (pre- 
sumably random) potential configurations, rather than 
calculating average quantities. We first focus on the 
almost depleted wire and recover the earlier results of 
Refs. [12l [17] for average phase boundary for the effec- 
tive p-wave model. We particularly recover the result 
that the disorder is always detrimental to the topolog- 
ical order for p-wave superconductors. We then show 
how for individual disorder configurations, one can re- 
late the phase diagram to an experimentally accessible 
quantity: the normal state conductance. This result al- 
lows us to solve inter alia the gaussian disordered p-wave 
problem exactly for all values of the disorder strength. 
Next, we focus on the (superconducting) semiconducting 
wire at arbitary doping and show that the topological 
phase transition at higher doping is it not described by 
the effective models. Most importantly, the disorder is 
not detrimental to topological order, rather the topologi- 
cal region is shifted to higher chemical potentials and its 
area in phase space is conserved. Strikingly, if the scat- 
tering is regular e.g. due to a superlattice, the area of the 
topological phase can even be made to increase beyond 
the clean value. 

While our main aim is to study the topological phase 
transitions in the semiconducting wire/s-wave supercon- 
ductor system, we first focus on the spinless p-wave 
Hamiltonian. In the latter case, the calculation is eas- 
ier to follow, but still illustrates the essential concepts. 
Moreover, the p-wave model arises as the effective Hamil- 
tonian in the large B limit of the s-wave Hamiltonian. 
Comparing the p-wave with the full s-wave calculation 
will then show where the effective Hamiltonian approxi- 
mation breaks down. We note that the p-wave model was 
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solved for chemical potential set to the band center and 
for specific position-dependent potentials [18l[T9]. Here, 
we will present a general solution. 

The Bogoliubov-de Gennes (BdG) Hamiltonian of a 
spinless p-wave superconductor in Id is given by: 

H = h{p,x)Tz ^upTx, (1) 
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where h{p^ ^) ~ 2rn ^ V{x) — jj. is the (spinless) single 
particle Hamiltonian, p is the momentum operator, m 
the electron mass, V{x) an arbitrary scalar potential, /i 
the chemical potential, and u = A/pp with A the super- 
conducting gap and pp the Fermi momentum. Here and 
below Ti {i = X, z) denote the Pauli matrices in the 
electron-hole space. In order to make use of the chiral 
symmetry of the Hamiltonian, we first apply a global ro- 
tation in the electron- hole space (r^ ^ r^^, Tx Ty) and 
cast the Hamiltonian into off-diagonal form [20 . The 
main use of this form is that it is now easy to see that the 
zero mode solutions, i.e. Majorana fermion solutions, are 
either of the form x+ = ( "^q^ ) or of the form X- — )•> 
with {h{p^ x) ± iup) (f± = 0. The linear in momentum 
term can be removed by a gauge transformation with 
a suitably chosen imaginary parameter (p± = e^^'^^ip, 
where ku = mu/h. We then find that ip satisfies 

We identify this equation as the normal state equation 

with an effective chemical potential fi = ji with 

one crucial distinction: it is e^^^^^p that needs to be 
normalized, rather than ip itself. Thus diverging solutions 
ofEq. ^ as X ^ ±00 lead to normalizable wavefunctions 
^±7 provided the divergence is not faster than e^^^^. 

For the sake of concreteness we focus on an half infi- 
nite {x > 0) wire, i.e. we assume that at points x < is 
the vacuum state (a normal insulator), specified by the 
boundary condition x|o = (it is easy to generalize to 
boundary conditions of the form ax(^o)+^^Uo = 0) 
X is normalizable, i.e. X ^ sufficiently fast as x ^ 00. 
From standard Sturm-Liouville theory, recall that if the 
solutions of the (spinless) Hamiltonian ([2| are localized, 
then there is one exponentially decaying solution (which 
we choose to be /) and one exponentially increasing solu- 
tion (which we choose to be g) for large x. If the spinless 
electron is delocalized then both / and g are oscillatory. 
We choose a suitable linear combination ip = Af -\- Bg 
such that '^(O) =0 and hence also x fulfils the boundary 
condition. Then for large ip ^ e^^ with A real ^22^ and 
a function of the effective chemical potential A = A(/i). 
We identify three cases (i) A < — /c^, (ii) |A| < ku, and 
(iii) ku < A. For case (i) there are two zero modes x+ 
and X- =• This can only happen if the decaying solu- 
tion / itself accidentally fulfills the boundary condition, 
and the two solutions will be lifted away from zero for 




Figure 1: Topological charge Q = det(r) of a disordered p- 
wave nanowire as a function of chemical potential /x and dis- 
order strength 7, for a single disorder configuration in a short 
wire (L = 100a, with a the lattice constant). The inset shows 
a single disorder configuration in a long wire (L = 10000a). 
The red solid line in the main plot is the phase boundary 
computed from Eqs. (|3| and the normal state conductance C, 
the red solid line in the inset/red dashed line in the main plot 
from Eqs. ([3| and Q. The numerical calculation was done in 
a TB model with ku — 10a~^ and a chemcial potential in the 
leads /xiead = 0.5h^ /2ma^ . 

small perturbations, i.e. are not topologically protected. 
This case corresponds to an accidental crossing of energy 
levels at zero energy p3]. In case (ii) there is only one 
Majorana state, X- which is the topologically protected 
state, and in case (iii) there are no zero modes and thus 
no topological state. We thus obtain a formula for the 
topological charge: 

Q = sgn(n|A(/i - mu^/2) \ - mu). (3) 

This the central result for the p— wave part of our work. 

We are now at a position to demonstrate the topolog- 
ical robustness of the zero energy solutions. First note 
that it is only the asymptotic limit of the solutions tp of 
the effective Schrodinger equation that matters for the 
existence of the solutions. Next notice that local pertur- 
bations of the potential (unless infinite) cannot change 
the asymptotic limit of the solutions regardless of their 
size and shape. Thus if there is a zero mode of the 
BdG hamiltonian for some potential profile (i.e. it is in 
the topological state) so will any other Hamiltonian that 
differs from the former by a local perturbation, demon- 
strating topological invariance. 

For a disordered (normal-state) wire, A is usually 
called the Lyapunov exponent and can be estimated from 
the conductance as: A = — (2/L) log(G/Go), where L is 
the wire length and Go the conductance quantum [24 . 
Hence, for fixed Eq. (|3| allows one to determine the 
topological charge of a p-wave quantum wire from its nor- 
mal state conductance alone. In short wires A fluctuates 
strongly as the chemical potential varies, and as a conse- 
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quence there are multiple changes of the topological prop- 
erties. This is shown on the example of a single disorder 
realization in a short wire in Fig. [l] where we computed 
the topological charge within a tight-binding (TB) model 
numerically from Q = det(r) where r is the reflection ma- 
trix [25] . The topological phase boundary computed from 
Eq. (|3| and the numerically computed normal state con- 
ductance agrees very well with the det(r)-criterion; small 
deviations of the exact position of the phase boundary are 
due to finite size effects. 

For longer wires the Lyapunov exponent is a self aver- 
aging quantity, i.e. A(L) ^ A, as L ^ oo, where A is the 
average Lyapunov exponent. For a wire with gaussian 
disorder {V{x)V{y)) = 7^(0:: — y) at energy e, it can be 
obtained in closed form [26l [27] : 



A(e) 



(4a) 



l dln(Ai(-2V3x)^ + Bi(-2V^a:)^) 
= "2 Tx • 

Then the topological transition condition Eq. (|3| be- 
comes fi|A(/i — mu'^/2)\ = mu, valid for the entire range 
of /i, li, 7 and shown as a red solid line in the inset of 
Fig. [ij The inset also shows numerics for a single dis- 
order configuration for a long wire, demonstrating that 
due to the self- averaging long wires have a well-defined 
universal topological phase (similar numerics, but aver- 
aged over disorder was shown in [28 ). At high ener- 
gies, we have the golden rule result A ~ l/4^tr, where 
£tr = {iJL — mu^ / 2) / is the transport mean free path. 
We then obtain the condition that there is a topological 
transition at ku^-tr = 1/4, in agreement with Ref. [12 , 29 . 

From Eq. ([3| it can be also concluded that for /i > 
any scattering is detrimental to the topological phase: 
Then A = in the clean case (the normal state solutions 
are extended), and any scattering leads to A > 0. For Jl < 
topology can be in principle induced as seen in the inset 
of Fig. [ij There, a topological phase is created for /i < 
and 7 > due to states in the Lifshitz tail below the 
band bottom. This however is a relatively small effect. 
We shall see below this picture is drastically different for 
the experimentally relevant proximity nanowire systems. 

We now focus on the experimentally more relevant sys- 
tem: a nanowire with Rashba spin-orbit coupling in prox- 
imity to an s-wave superconducter. The BdG Hamilto- 
nian is then given as ^ |4j : 



H = h{p, x)tz + apayTz + Bax + Ar^,, 



(5) 



where h{p^ x) = /2m + V{x) — /i is the (spinless) single 
particle Hamiltonian, a the spin-orbit coupling strength, 
B the Zeeman splitting and A the induced s-wave order 
parameter. {i = x^y^z) are the Pauli matrices in spin 
space. The topological state appears for 5^ > A^+/i^. In 
this single orbital mode limit, the system is in class BDI, 



which is distinguished from class D by the presence of the 
chiral symmetry. This allows to bring the Hamiltonian 
into off-diagonal form [30] , and the zero-energy Majorana 
states are of again of the form x+ = ( ) X- = ( ^ ) • 
After a rotation by 7r/2 in a space around the x-axis and 
premultiplying with icr^^, we find that if± satisfies a 2 x 2 
nonhermitian eigenvalue problem: 



(h{p, x)cFz ±B ± Aax - iapax)(p± = 



(6) 



We now construct the zero energy solution for small a. 
Larger values do not change the picture, but rather renor- 
malize the topological-normal phase boundaries. We first 
let (p e~'^^(p^ where is an order a parameter that is 
yet to be determined. Then we have p ^ p ^ ifin. Next 
we collect terms of order a and treat them as perturba- 
tions. We then have H = Hq + i^i. 



Ho = h{p, x)crz + 5 + Ada, 
. hnp 



(7a) 



Hi = —iapcFx + i '-—^(Jz + hnacFx cr^ • (7b) 

m 2m 

The last two terms can be absorbed into Hq by redefin- 
ing jl and A. The eigenfunctions of Hq are of the form 
^±'0(^; ^) where is the local solution of h^l) = ei/j^ and 
^±{e) is the eigenspinor of the 2x2 matrix ea^ + Aa^^. 
We now choose hz such that Hi anticommutes with this 
matrix. Then it is easy to see that 

+e+(-e)e'^^(C/(x; -e) + Dg{x; -e)), (8) 



where e = — A^ and k, = maA/he^ is a local solu- 
tion of {Hq -\- Hi)(f = to order a^. As in the p-wave 
case we have written as a sum of the two linearly inde- 
pendent solutions / (decaying) and g (increasing). Then, 
(p-^ is a valid zero-energy solution (and thus a Majorana 
fermion) if it is normalizable and satisfies the boundary 
conditions. Repeating the calculation for the the other 
sector {(p-)^ we obtain n ^ —n. 

We assume again without loss of generality that the 
system is in a normal insulator state for x < and 
the boundary condition (p{{)) = (note that is a 
spinor). We then identify three cases: (i) If 5 > A, and 
|A(/i±e)| < \f<i\ or |A(/i±e)| > there are two decaying 
and two diverging solutions. Then the boundary condi- 
tion at X = cannot be satisfied except in accidental 
cases such that /(±e) already fulfil the boundary condi- 
tion. Then there is also a second solution in the other 
sector, and the zero-energy states are not protected. The 
system is thus in the trivial state with the possibility of 
accidental zero modes, (ii) If 5 < A, then both k and e 
are imaginary, hence there are always two decaying and 
two diverging solutions. However, there are no acciden- 
tal zero modes with /(±e) already fulfilling the boundary 
condition because this would mean / is an eigenfunction 
of (Hermitean) h with an imaginary eigenvalue, (iii) If 



B > A and |A(/i ± e)| < \k\ < |A(/i T there are 
one diverging and three decaying solutions in one sector 
and one decaying and three diverging solutions in the 
other sector. Then the boundary condition at x = can 
be generally satisfied in the sector that has three decay- 
ing solutions and there is a Major ana state. As before, 
the solution is robust, because local perturbations do not 
change the asymptotic behavior of / and g. In summary 
we have: 

^ A A/ M maA\ /, ^ , maA\ 

Q = sgn (^|A(/i + e)| - j sgn (^|A(/i - e)| - j 

(9) 

This is our central formula for the s-wave case. Note 
that the first term in Eq. ([9| reduces to Eq. (|3| in the 
large 5-limit; in this case the second term induces new 
physics. 

We now discuss various types of scattering using 
Eq. ([9| and start with the example of a superlattice. In 
the clean case, the required odd number of channels for 
the topological state is only achieved if the chemical po- 
tential is within the Zeeman gap and hence close to the 
band bottom [3[ H]. The superlattice allows this for a 
larger range of chemical potentials out of the Zeeman gap. 
The minigaps (or equivantly, perfect backscattering) can 
induce topological order away from the band bottom, and 
the topological phase space area can be even increased 
as compared to the clean case, as shown in Fig. [2](a,b). 
Strikingly, one can even use topologically trivial pieces 
to induce the topological state. 

In the experimentally relevant case of irregular scat- 
tering, we use the average Lyapunov exponent given by 
Eq. Q to determine the phase boundaries of a long quan- 
tum wire from Eq. ([9|. Noting that A is a monotonous 
function of energy, we get: 

/i± = F-\m'^'XaA/^/B^ - A^)/\^ ± ^ - A^ (10) 

In the weak disorder limit, A ^ oo, we recover the clean 
wire result :/i± = — A^. In contrast to the com- 

mon wisdom based on the effective p-wave model, we find 
that the topological region is not destroyed by disorder 
but merely shifted to higher chemical potentials. In fact 
the chemical potential (or gate) range where the wire is 
topological, — /i_ = 2\/B'^ — A^, is independent oi the 
disorder strength, while the total area of the topological 
region in the (5,/i) plane is conserved. We stress that 
this result is valid to all orders in disorder strength. 

Fig. [2|c,d) shows our numerical results for a nanowire 
with disorder. In agreement with our analytical results, 
we observe that for a long wire the topological phase is 
merely shifted to larger chemical potentials. The disor- 
der then creates a well-defined topological region for a pa- 
rameter range where the clean wire is trivial. In fact, dis- 
order is even favorable for experiments to form Major ana 
fermions, as the threshold magnetic field is considerably 
lowered for ji away from the band edge. In a short wire. 
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Figure 2: Topological charge Q = det(r) as a function of chem- 
ical potential ji and Zeeman splitting B for a (a) clean system, 
(b) a superlattice, and (c, d) disorder. Red lines in (b-d) are 
phase boundaries calculated from Eq. green dashed lines 
show the clean phase boundary for comparison, (b) The su- 
perlattice (see inset) parameters were 

A — T? I2m\? ^ and fco — ^.^^h~^ with fco — ma/h^ and 
the numerical calculation was done using a transfer matrix 
method in Mathematica. The numerical calculations in (c- 
f) were done within a TB model: (c) shows the topological 
charge for a single disorder realisation in a short (L = 100a 
with a the lattice constant) and (d) in a long (L = 4000a) 
wire, (e) and (f) the respective tunnel conductances for a 
fixed II — 0.3t, with t — fi? /2mo?. White dashed lines in 
(e, f) indicate the boundaries of the topological phase in 
(c, d). The remaining TB parameters were /cso = 0.05a~^, 
A = 0.15t, 7 = 0.06t^ , and the chemical potential in the 
leads /iieads = 0.5t. For the tunneling conductance in (e, f) 
a barrier of height 1.5t was added on one lattice site next to 
one end of the wire. 



the topological phase is more fragmented due to the fiuc- 
tuations in the normal state conductance in agreement 
with Eq. ([9|. Nevertheless, both for a short and a long 
wire, a clear Majorana zero-bias peak (ZBP) appears in 
the tunneling conductance, as shown in Fig. [2|e,f). We 
note that the wire would be in the trivial phase without 
Majorana fermions without disorder for the whole range 
of parameters shown in Fig. [2|e,f). 
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Recently, it was argued that ZBPs in nanowires may 
appear even without Majorana fermions [3TH33] . Here 
we caution against this interpretation. As a ZBP out of 
the clean topological phase boundary may well be a Ma- 
jorana fermion within the dirty topological phase bound- 
ary, in particular if B > A and the ZBP remains for a 
range of magnetic field. In fact, for a dirty wire all ac- 
cidental zero mode solutions will shift under changes in 
the magnetic field. 

In conclusion, we studied the effects of single-particle 
scattering on the topological properties of a quantum 
wire in contact with an s-wave superconductor. We ob- 
tained analytical formulas for the phase boundaries valid 
for regular and irregular scattering. In the case of ir- 
regular scattering, our formulas apply to all orders in 
the disorder strength as well as to individual members 
of the ensemble. Our main result is that disorder does 
not destroy topological order, in contrast to the general 
expectation based on effective models, rather the topo- 
logical phase is merely shifted to higher chemical poten- 
tials, while the total phase area at fixed magnetic field 
is conserved. Moreover one can even increase the topo- 
logical phase area with periodic modulation of the gate 
potential. 
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